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We perform no-core full configuration calculations for the Lithium isotopes, 6 Li, 7 Li, and 8 Li 
with the realistic nucleon-nucleon interaction JISP16. We obtain a set of observables, such as 
spectra, radii, multipole moments, transition probabilities, etc., and compare with experiment where 
available. We also present one-body density distributions for selected states. Convergence properties 
of these density distributions shed light on the convergence properties of one-body observables. We 
obtain underbinding by 0.5 MeV, 0.7 MeV, and 1.0 MeV for 6 Li, 7 Li, 8 Li respectively. Magnetic 
moments are well-converged and agree with experiment to within 20%. 
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I. INTRODUCTION AND MOTIVATION 

The rapid development of ab initio methods for solv- 
ing Unite nuclei has opened a range of nuclear phenom- 
ena that can be evaluated to high precision using realis- 
tic inter-nucleon interactions. Here we investigate three 
Li isotopes through direct solution of the nuclear many- 
body problem with the JISP16 nucleon-nucleon (NN) in- 
teraction [H-III by diagonalization in a sufficiently large 
basis space that converged binding energies are accessed 
by simple extrapolation. The many-body approach we 
adopt is referred to as the no-core full configuration 
(NCFC) method which yields uncertainty estimates 
for binding energies. 

We adopt the traditional harmonic oscillator (HO) ba- 
sis which enables us to isolate and remove spurious effects 
of center of mass (cm) motion from all observables and 
from the one-body density matrices (OBDMs) to high 
precision. Since we address light nuclei, we feel this ca- 
pability is an important ingredient. A further advantage 
in using the (HO) basis is its ease in performing analytical 
evaluations and straightforward matrix element calcula- 
tions though certain alternative basis choices also have 
these advantages @, I n this context, the NCFC ap- 
proach is similar to the no-core shell model (NCSM) J9|. 
The main differences are that in the NCFC approach we 
do not use the Lee-Suzuki renormalization procedure [l(| 
which is commonly employed in the NCSM; and more 
importantly, we estimate the numerical accuracy of our 
results based on the rate of convergence and dependence 
on the basis space parameters 0, Q . 

We use a set of finite single-particle HO bases, charac- 
terized by two basis space parameters, the HO energy Ml 
and the many-body basis space cutoff N max . -/V max is de- 
fined as the maximum number of total oscillator quanta 
allowed in the many-body basis space above the min- 
imum for that nucleus. Independence of both param- 
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eters hQ and ./V max signals numerical convergence; for 
bound states, true convergence is generally expected to 
be reached in the limit of a complete (infinite dimen- 
sional) basis. For the binding energy we use an extrap- 
olation to the complete basis space, with error bars re- 
flecting the uncertainty due to the extrapolation. 

In this work, we look at the actual nucleon densities, in 
addition to observables such as spectra, radii, and mul- 
tipole moments. In order to do so, we introduce a tech- 
nique for unfolding the cm motion from the one-body 
density. This allows us to obtain the translationally- 
invariant densities, without any smearing effects from the 
cm motion. Indeed, salient details of the density are of- 
ten enhanced in the translationally-invariant (ti) density, 
compared to the single-particle densities that are com- 
monly used in configuration interaction calculations, as 
we will show. 



II. METHODS AND INGREDIENTS 
A. Nuclear Hamiltonian 

We begin with the translationally-invariant Hamilto- 
nian for the A-body system in relative coordinates 



v = 1 (Pi-Pj) 2 

A 2m 

t<j 



(1) 



l< 3 



where mc 2 is the nucleon mass times the speed of light (c) 
squared (taken to be 938.92 MeV, the average of the pro- 
ton and neutron masses), Vmn is the NN interaction, and 
Vc is the Coulomb interaction, which acts between the 
protons only. We adopted the NN interaction JISP16, a 
realistic NN interaction initially developed from NN data 
using inverse scattering techniques. It is then adjusted 
with phase-shift equivalent unitary transformations to 
describe light nuclei without explicit three-body inter- 
actions [il-Q- 

JISP16 provides good convergence rates for the ground 
state (gs) energies of nuclei with A < 16. We investigate 
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convergence rates for a selection of additional observables 
in the present work, in particular the spectra, radii, and 
multipolc moments. We use the naive point-like oper- 
ators for these observables; we do not take effects such 
as meson-exchange currents into account. However, it is 
known that these effects can cause significant corrections 
to observables such as magnetic moments . We should 



therefore expect similar deviations from experiments for 
our calculations of these observables. 



B. No-Core Full Configuration calculations 

In the many-body framework that we are using, we 
expand the nuclear wavefunction $ in a basis of Slater 
determinants of single-particle HO states. Note that we 
use single-particle coordinates, rather than relative coor- 
dinates, in the nuclear wavefunction. That means that 
our wavefunctions, and therefore our OBDMs calculated 
as expectation values of one-body operators, will include 
cm motion. 

All many-body basis states are included with HO 
quanta up to and including the amount governed by 
the iV max truncation. Thus, if the highest HO single- 
particle state for the minimal HO configuration has Nq 
HO quanta, then the highest allowed single-particle state 
in the truncated basis will have Nq + N max HO quanta. 
Furthermore, our calculations are 'No-Core' configura- 
tion interaction calculations. This means that all nuclc- 
ons participate in the interactions on an equal footing. As 
we increase A max , and approach convergence, we expect 
physical observables to become independent of both the 
HO parameter HQ and the truncation parameter A max . 
However, due to current limits to our finite basis, our 
calculations do show some parameter dependence, even 
in the largest basis spaces. We apply previously estab- 
lished extrapolation tools to take the continuum limit of 
the binding energy as we discuss shortly. 

The HO basis for single-particle states, in combina- 
tion with this many-particle A max truncation, leads to 
exact factorization of the nuclear wavefunctions into a 
cm wavefunction and a ti wavefunction 



*(rj) = $? m (i?) ® </> ti 



(2) 



where R = (4;) 2j=i anc ^ 4>ti depends only on inter- 
nucleon coordinates. The actual Hamiltonian that we are 
using involves adding and subtracting H cm so that Ha 
takes the form 
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In order to separate the cm excited states from the low- 
lying states of interest, we adopt the Lawson method [l2| 
whereby we add a Lagrange multiplier term, X(H^ m — 



§7i.f2), to the many-body Hamiltonian, Eq. ([3]) 
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With A positive, states with cm excitations are sep- 
arated by multiples of Xhtt from the states with the 
lowest HO cm motion. Since the Lagrange multiplier 
term acts only on the cm coordinate, it is independent of 
the inter-nucleon coordinates and it does not affect the 
energy eigenvalues or the translationally-invariant wave- 
functions 4>ti of the low-lying states. Indeed, observables 
for the low-lying states are independent of A, as long 
as Xhfl is much larger than the excitation energy of the 
highest state of interest. 

In the truncated basis space, we can now write the 
many-body Schrodinger equation as a finite matrix equa- 
tion with a real, symmetric, sparse matrix. The eigen- 
values of this matrix give us the binding energy, and the 
corresponding eigenvectors give us the wavefunctions. In 
any finite basis space, the eigen-energies satisfy the vari- 
ational principle and show uniform and monotonic con- 
vergence from above with increasing A max , allowing for 
extrapolation to the infinite basis space. To obtain the 
extrapolated gs energy E gs (oo), we use a fitting function 
of the form 

-B 5S (A max ) = a exp(-c N max ) + E gs (oo) . (5) 



This is an empirical method [4]46| that is valid within 
estimated uncertainties that we now define. We assign 
equal weight to each of three successive values of A max 
at a fixed HQ and perform a regression analysis. The 
difference between extrapolated results from two consec- 
utive sets of three A max values is used as the estimate of 
numerical uncertainty associated with the extrapolation. 
The optimal HQ value for this extrapolation appears to 
be the hn that minimizes the difference between the ex- 
trapolated energy and the result at the largest A max . 
Typically, this corresponds to a hQ value slightly above 
the variational minimum. Of course, the extrapolated 
results should be independent of hQ, within their numer- 
ical error estimates, and we do check for such consistency. 
Furthermore, we often adjust our numerical error esti- 
mate by considering the results over a range of 5 MeV in 

nsi. 

For other observables, we do not have a robust and re- 
liable extrapolation method; we therefore use the degree 
of (in) dependence from the basis space parameters Ml 
and A max as a measure for convergence as we describe 
further below on a case-by-case basis. 

For nuclei with mass A > 4, it is challenging to ob- 
tain convergence (independence of both A max and f2) for 
all observables within the NCFC approach with realistic 
interactions. As N max increases, the dimension of the 
many-body basis increases exponentially, and there is a 
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clear need for high-performance computing. For the Li 
isotopes we investigate here, we have limited the basis 
space to about 1 billion: for 6 Li the largest basis space 
is N max = 16, with a dimension of 805,583,856; for 7 Li 
the largest basis space is N max = 14, with a dimension 
of 1,244,131,981; and for 8 Li the largest basis space is 
iV max = 12, with a dimension of 1,222,330,036. We use 
the code MFDn (Many Fermion Dynamics for nuclear 
structure) [H, [l4[ , which is a state-of-the-art numerical 
code for no-core configuration interaction calculations. 
The calculations were performed on the Cray XT4 and 
XE6 at the National Energy Research Supercomputer 
Center (NERSC), using up to 8,000 processors (4 or 6 
cores/processor) with 8 GB of memory each. 



C. One-Body Density Matrix 

The OBDM represents, in a compact form, sufficient 
information about quantum states of a system to eval- 
uate all observables that can be expressed by one-body 
operators. For example, a single wavefunction of 6 Li at 
iV max = 16 has a size of 5 GB, while the OBDM for a 
single state in that same basis is less than 5 MB. The 
OBDMs also provide the necessary information for visu- 
alizing nuclear density distributions as we demonstrate 
below. 

The non-local one-body density in coordinate space 
for an initial A-body wavefunction and final A-body 
wavefunction ^ f is defined as 

p fi {rx,f'x) = j mri,r2,...,r A ) 

x H t {r l ,r 2 ,...,r A ) d 3 r 2 ...d 3 r A . (6) 

The limit f\ = of Eq. (j6|) gives the local one-body 
density. For = ^ this corresponds to the probability 
of finding a nucleon at position r\ when the system is in 
that state. That is, the local one-body density distribu- 
tion is given by 

pS(n) = [p(nrt)S(n-i*i)d 3 r[ (7) 



where we suppress the state labels for simplicity and in- 
sert a superscript il to signify the dependence on the HO 
basis space used in the evaluation of the eigenfunctions. 
Since it depends on the single-particle coordinates, we 
refer to this density as the space-fixed (sf) density. 

Note that, due to our use of single-particle coordinates, 
rather than relative coordinates, our wavefunctions ^f(fi) 
include cm motion. The resulting one-body density dis- 
tributions will therefore include contributions from the 
cm motion. However, because of the exact factoriza- 
tion of the cm wavefunction and the ti wavefunction (see 
Eq. ([2])) this density is actually a convolution of the cm 
density p^ m and the ti density Pti(r) describing the prob- 
ability of finding a nucleon at position f relative to the 



cm of the entire nucleus 1 

Psdn) = I Pti(n - R) p'L(R) d 3 R. (8) 



For the HO basis, p^ m is a simple Gaussian (the gs density 
of H cm ) with explicit dependence on Q that smears out 
pti- This smearing can obfuscate interesting details of pu- 
Furthermore, it introduces a spurious dependence on the 
basis parameter il into p^ f that masks the convergence. 
Even in the limit of a completely converged calculation, 
the single-particle density p^ { depends on O, whereas p t i 
becomes independent of the basis. 

In order to eliminate these smearing effects and to help 
develop a physical intuition for the ah initio structure of 
a nucleus, it would be helpful to see the coordinate space 
density distributions free of spurious cm motion. This 
can be achieved by a deconvolution of the cm density 
and the ti density using standard Fourier methods [l5| 



Pti(ri) 



F~ 



F[p%(n)} 
F[ P ? m (R)] 



(9) 



where F[f(r)] is the 3-dimensional Fourier transform of 
f(r). At convergence, the dependence on f2 should cancel 
on the RHS of this equation. That means that after this 
deconvolution, we can better investigate the convergence 
of the physically-relevant ti density. 

In addition to the 3-dimensional ti densities, we also 
show results for different multipoles of the ti density, 
p[^\r), defined with the same initial and final states, 
having total angular momentum J and magnetic projec- 
tion M, as 



2.1 



K=0 



(JMK0\JM) 
V2J + 1 



YK{?)p\?\r). (10) 



These multipoles allow for a better assessment of the nu- 
merical convergence of the densities, as we will illustrate 
below. Another advantage of the multipolc expansion is 
that it allows for a straightforward determination of the 
ti density for any M: the multipoles p[^\r) arc inde- 
pendent of M. However, some features of the density, 
in particular clustering, are more apparent from the 3- 
dimensional plots of p{r) than from plots of the different 
multipolc components. 

The results we present here for densities of the Li- 
isotopes are all translationally-invariant unless we state 
otherwise. An alternative method to isolate the ti density 
has been developed in [TtJ ■ 



1 Note that our ti density is not the same as the one-body den- 
sity in Jacobi coordinates, cr(£) as defined in Ref. [H, which 
describes the probability of finding one nucleon at £ relative to 
the cm of the remaining A — 1 nucleons. Both these densities are 
translationally-invariant, and they are related to each other via 
rcscaling factors 
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D. Observables 

In HO space, the space-fixed OBDM is specified by its 
matrix elements 



fi 



(ii) 



where a and (3 stand for a set of single-particle quantum 
numbers (n a , l a , j a , m a , r„, a ) and (np,lp,jp,mp,T z> p), 
and we use the Dirac bra-ket notation to represent the to- 
tal many-body state vector. Once we have obtained the 
one-body density matrix elements pp a (OBDMEs), we 
can easily calculate observables that can be expressed as 
one-body operators. For initial and final states with to- 
tal angular momentum Jj j and possibly additional quan- 
tum numbers Xij, but with the same magnetic projection 
M, the E2 matrix elements using the canonical one-body 
electromagnetic current operator are given by 



Mi 



fi 
E2 



(\fJfM\E2\XiJiM) 

E 4a M / r 2 Y 2 °(f)d 3 r[ 

a/3 



(12) 



with a and j3 restricted to the protons only (t z = ^). 
Note that the fact that the OBDM includes cm motion 
does not matter for E2 matrix elements (nor for Ml ma- 
trix elements discussed below): the cm wavefunction is a 
normalized s-wave, and does not contribute to the inte- 
gral due to the factor Y^f). 

For comparison with experiments, it is more conve- 
nient to convert these M-dependcnt matrix elements to 
reduced matrix elements using the Wigner-Eckart theo- 
rem [l8|. For a proper tensor operator T^k the reduced 
matrix element is defined by 



(XfJfWTKWXiJi) = (XfJfMflTKklXiJiMi, 



(J f M f Kk\JiMi] 



(13) 



provided that the Clebsch-Gordan coefficient in the de- 
nominator, (JfMfKk\JiMi) (following the conventions 
of Ref. Pil), is not zero. In terms of the reduced E2 
matrix elements, reduced E2 transition probabilities are 
given by [l9| 



B(E2;»^/) 



2J, 



— -(XfJf\\ E2 p^) 2 (14) 



in units e 2 fm . The quadrupole moment is convention- 
ally defined through the E2 matrix element for M = J 



16tt 



1/2 



Q = y—J ( XJM = J \ E2 \ XJM = J ) (I 5 ) 

and can also be expressed in terms of the reduced matrix 
element as Oil 



1/2 (JJ2QL7J) 
V2J + 1 



(AJ||E2||AJ) (16) 



in units e fm 2 . 

The matrix elements for the Ml transitions and mag- 
netic moments receive contributions both from the pro- 
ton and neutron intrinsic spins and from the proton or- 
bital motion. Again, we consider only the canonical 
one-body electromagnetic current operator, in which case 
they can be calculated from the OBDMEs 



M fi 

IU M1 



(XfJ f M\Ml\XiJiM) 
= E 4° (a\kl + T z )(L + g P a) + \{l 



where g p = 5.586 and g n = —3.826 are the proton and 
neutron gyromagnetic ratios in nuclear magneton (/j-n) 
units; the quantities L, a and r represent the conven- 
tional orbital angular momentum, spin and isospin oper- 
ators. In terms of the reduced Ml matrix element, the 
reduced Ml transition probabilities are given as [19( 



B(Ml;i->/) = 



1 



2.L 



j(X f J f \\ Ml ||A 4 J ? ) 2 (17) 



in units /ify, and the magnetic moment is defined as 

'4^ 1/2 (^10|JJ) (AJ||M1||AJ) (18) 



1 3 J V2JTT 



in units /ijy. 

We also present results for the point-proton root-mean- 
squarc (RMS) radius, (r^) 1 / 2 . This can be calculated ei- 
ther from the translationally- invariant local density, pti, 
or directly, as a two-body operator, from the wavefunc- 
tion in single-particle coordinates, ^(fl). We performed 
the calculations with both approaches and confirmed the 
results were identical to within the numerical precision 
of four significant figures. 

In order to convert measured nuclear charge radii, 
(r 2 ) 1 / 2 , to point-proton radii, we use [13] 



c>2 N 2 

— a„ —tt~ — 



3r 



4M 2 c 2 



(19) 



Here, i? 2 = 0.769(12) fm 2 is the RMS proton charge ra- 
dius, i? 2 = -0.177(4) fm 2 the RMS neutron charge ra- 
dius, M p the proton mass, and 4 ^ 2c 2 ~ 0.033 fm 2 the 
Darwin-Foldy correction. 

In addition to these correction terms there is also a 
spin-orbit contribution to the charge radius, but this 
contribution is model-dependent, and (for the nuclei dis- 



cussed here) less than 1% for realistic wavefunctions [21 
Hence we neglect it. For the experimental radii we 
use the values of Ref. !2j3], which were obtained from 
high-precision laser spectroscopy measurements of iso- 
tope shifts in combination with the 6 Li charge radius as 
absolute reference. Using the 7 Li charge radius as ref- 
erence in combination with the same isotope shifts gives 
radii that are about 2% to 3% smaller than the ones we 
have listed in the tables below [221. 
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FIG. 1: (Color online) The gs energy (top) and excitation 
spectrum (bottom) of 6 Li for a sequence of JV max values (in- 
dicated in the legend) and as function of the HO energy. The 
extrapolated gs energy is shown at specific values of HQ with 
uncertainties (defined in the text) indicated as error bars. 



III. RESULTS 
A. 6 Li 

In Fig. [TJ we compare the gs energy and excitation 
energies at a sequence of iV max values and as a function 
of the HO energy hil. We also provide the extrapolated 
gs energy as a function of hil along with the assessed 
uncertainties (error bars) as described above. 

The gs energy for 6 Li is rapidly converging as indicated 
by the emerging independence of the two basis parame- 
ters (-/V max , hfl). The convergence is most rapid around 
fiil = 17.5 to 20 MeV, where the variational upper bound 
on the energy is minimal. Our extrapolated gs energy [B| 
shows that the system is underbound by 0.50 MeV. Ex- 
citation energies are well converged at higher iV max (12 
and above) values, at least for 3 + and + states. Note 
that these states are narrow resonances: the experimcn- 
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FIG. 2: (Color online) The RMS point-proton radius of the 
gs of 6 Li as a function of HO energy at various iV max values 
(top) and as function of ./V max at various values of the HO 
energy (bottom). 



tal width of the 3 + is 24 keV and the width of the first 
excited + is 8 eV. 

The excitation energy of the first 2 + state is much less 
converged, and shows a systematic increase with increas- 
ing HQ. Such ^-dependence of the excitation energy is 
typical for wide resonances as observed in comparisons 
of NCSM results with inverse scattering analysis of a- 
nucleon scattering states [HHH . In light of these previous 
analyses, the significant ?if2-dependence seems commen- 
surate with the large experimental width of 1.3 MeV for 
this 2+ state. 

In Fig. [21 we show the dependence of the RMS point- 
proton radius on the basis space parameters iV max and 
hfl for 6 Li. It appears that this radius is converging 
less rapidly than the gs and excitation energies. Further- 
more, the convergence is neither monotonic nor uniform: 
at small values of the HO energy the radius tends to de- 
crease with increasing -/V max , whereas at larger values of 
Hfl the radius increases with increasing -/V max ; around 
hil = 10 to 12.5 MeV the RMS radius is nearly indepen- 
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FIG. 3: (Color online) The angle-averaged density of the Li 
gs for various iV max values at hfl = 10 MeV (blue curves) 
and 17.5 MeV (red curves) on on a linear (top) and semi- 
logarithmic (bottom) scale. 

dent of N max . Because of this, it is difficult to make re- 
alistic estimates of uncertainties for radii and other long- 
range obscrvablcs such as quadrupole moments which ex- 
hibit similar patterns. 

The convergence patterns shown in Fig. [2] may be un- 
derstood from the following observation: since the HO 
wavefunctions fall off like a Gaussian, e~ cr while the true 
nuclear wavefunction falls off like an exponential, c~ dr , 
observables whose calculations are weighted towards the 
tail of the wavefunction, such as the RMS radius, will 
converge slower than those obscrvablcs that depend less 
on the tails, such as the energy and the magnetic mo- 
ment (see below) . Furthermore, it is well known that the 
RMS radius, and also other long-range operators such as 
the quadrupole moment, are minimally affected by the 
short-range correlations 0, . And the value of hfl 
that minimizes the gs energy is not necessarily the value 
of frfl that best represents the long-range behavior of the 
wavefunction. 

In Fig.[3]we show the radial density distribution for two 



sets of finite basis spaces, with TiQ, = 10 and 17.5 MeV 
respectively. The lower panel shows that the exponen- 
tial tail is much better represented in a HO basis with 
hfl = 10 MeV than in a HO basis with 17.5 MeV, and 
that the long-range behavior of the one-body density is 
therefore converging much more rapidly in a HO basis 
with hQ = 10 MeV, and rather poorly converging in the 
HO basis that minimizes the gs energy. That is, the ra- 
dial density calculated with hil = 10 MeV shows much 
more consistent long-range behavior at the three highest 
N max values than the density calculated with hQ = 17.5 
MeV. This leads us to the conclusion that while the value 
of HQ that minimizes the gs energy is an appropriate 
value when calculating gs and excitation energies, as well 
as magnetic obscrvablcs, this value is not necessarily opti- 
mal for calculations of obscrvablcs that depend primarily 
on long-range correlations, even in the moderately large 
basis spaces considered here. 

Therefore we will quote results for long-range observ- 
ables at the h£l value where the RMS radii for various 
Amax values intersect as seen in the top portion of Fig. [51 
rather than at the h£l value that minimizes the gs en- 
ergy. To be specific, for the Li-isotopes under discussion 
here, we simply take the results at Ti£l = 10 to 12.5 MeV 
(where the N ma _ x dependence appears to be minimal) as 
our approximation to the converged value of the RMS 
radius. In a similar fashion, we will cite results for the 
region of minimal N max dependence for other observables 
that depend primarily on long-range correlations. Such 
observables include RMS radii, E2 moments, and B(E2) 
transitions. Robust extrapolations to the infinite basis 
space and reliable error estimates for these observables 
remain an open question. 

In Fig. 3] we show the radial quadrupole component of 
the 6 Li gs density distribution, p^\r), at frfl = 17.5 MeV 
for several -/V max values. From this figure it is evident that 
the quadrupole moment 

Q = J P (2 Hr)r*dr (20) 

receives both positive and negative contributions. For 
JV max = 4, the positive contributions are significantly 
larger than the small negative contribution around 2 fm, 
but as N max increases, this negative region becomes more 
pronounced. At A r max = 16 the quadrupole density is 
nearly converged inside about 3 fm. The correspond- 
ing quadrupole moment, Q = —0.075 e fm 2 , appears 
to approach convergence. The small positive region in- 
side 0.9 fm contributes less that 1% to this value, but 
about half of the contribution from the negative region 
between 0.9 fm and 3.4 fm is cancelled by the positive 
tail of the quadrupole distribution: their contributions 
are —0.167 e fm 2 and +0.092 e fm 2 respectively. Note 
that this positive tail is well outside the charge radius. 

Table U presents a capsule view of selected spectral and 
other obscrvablcs for 6 Li. In principle, we can extrapo- 
late not only the gs energy, but also the binding energies 
of the excitated states to the infinite basis space. We 
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FIG. 4: (Color online) The radial quadrupole density p£ (r) 
of the 6 Li gs for various iV max values at hQ — 17.5 MeV on a 
linear (top) and semi-logarithmic (bottom) scale. 



can then calculate the excitation energy as the differ- 
ence between the extrapolated binding energies, treating 
the numerical error estimates as independent. For the 
excitation energies that appear to be converged (in par- 
ticular the 3 + state in Fig.Q]), such a procedure leads to 
an overestimate of the numerical uncertainty in the exci- 
tation energy: part of the numerical error is common to 
the ground state and the excited state. For such states 
we reduce the obtained numerical error estimates based 
on the apparent convergence of the excitation energies in 
finite basis spaces. On the other hand, for states with ex- 
citation energies that show a significant htt dependence 
(such as the 2+ state in Fig. [I}, we increase our numerical 



6 Li 


Expt. 


JISP16 AV18/IL2 CD-B 


E h (l+,0) 


31.994 


31.49(3) 


32.0(1) 


29.07 


( r pp) 


2.45(5) 


2.3 


2.39(1) 


2.25 


£x(3+,0) 


2.186(2) 


2.56(2) 


2.2(2) 




£x(0+,l) 


3.56(1) 


3.68(6) 


3.4(2) 




£ x (2+,0) 


4.312(22) 


4.5(3) 


4.2(2) 




£x(2+,l) 


5.366(15) 




0.0{Z ) 




Q(l+,0) 


-0.082(2) 


-0.077(5) 


-0.32(6) 


-0.066 


Q(3 + ,0) 










M (l+,0) 


0.822 


0.839(2) 


0.800(1) 


0.843 


^(3+ 0) 




1.866(2) 






B(E2;(3+ 0)) 


10.7(8) 


6.1 


11.65(13) 




B(E2;(2+,0)) 


4.4(23) 


7.5 


8.66(47) 




B(M1;(0 + ,1)) 15.43(32) 


14.2(1) 


15.02(11) 




B(Ml;(2+,0)) 




< 0.001 


0.002(1) 




B(M1;(2+,1)) 


0.1 (3) 


0.05(1) 






Mgt 


2.170 


2.227(2) 


2.18(3) 





2.14 



0.080 



0.843 



TABLE I: Selected Li observables calculated up through 
JV m ax = 16. The energies are in MeV; the RMS point-proton 
radius is in fm; the quadrupole moments are in e fm 2 ; the 
magnetic moments are in /xjv; the reduced B(E2) transition 
probabilities are in e 2 fm 4 ; and the reduced B(M1) transi- 
tion probabilities are in fx%. All listed transitions are to the 
ground state. The energies are obtained from extrapolations 
to the infinite basis space, with error estimates as discussed in 
the text; the dipole observables as well as the gs quadrupole 
moment are converged within the quoted uncertainty; the 
other quadrupole observables and the RMS point-proton ra- 
dius are evaluated at h£l = 12.5 MeV. We used Ref. [21)1 for 
the experimental value of the RMS radius and Ref. [26| for 
GT matrix element; the other experimental values are from 
Refs. [13, [H]. AV18/IL2 results are from Refs. [ll], [HQ 
and include meson-exchange corrections for the dipole observ- 
ables; CD-Bonn and INOY results are from Ref. [32] calcu- 
lated at N max = 16 and hil=ll and 14 MeV respectively, with 
the INOY gs energy extrapolated to the infinite basis space. 



error estimate based on this hil dependence. 

In general, magnetic dipole observables tend to con- 
verge rapidly. Indeed the magnetic moments of the gs 
and first excited state, as well as the B(M1) transitions 
to the gs, are well converged, with, at N msix = 16, a 
residual dependence on the basis space parameters that 
is of the same order as the overall accuracy of our lowest 
eigenvalue. 

For the RMS point-proton radius, as well as the B(E2) 
transitions to the gs, we list our results from the largest 
basis space, N niax = 16, at HQ = 12.5 MeV, as discussed 
above. The electric quadrupole moment of the gs is in ex- 
cellent agreement with the experimental value - perhaps 
better than might be expected due to basis space limita- 
tions and long-range nature of the quadrupole moment 
operator. 



8 



Both the gs binding energy and the excitation ener- 
gies calculated with JISP16 compare favorably to those 
calculated with alternative realistic NN interactions, i.e. 
Argonne V18^CD-Bonn JNOY, and SRG evolved N3LO 
interactions [5^, HO, HI, HH . However, with the addition 
of approp riate three-body interactions, such as Illinois- 
2 [30l l3l| or chiral three-body forces [34T - [36| , one can 
obtain somewhat better agreement with data than with 
the NN-only interaction JISP16. 

The obtained RMS point-proton radius is similar to 
that obtained with CD-Bonn: both are about 10% too 
small compared to experiment. INOY [32[ gives an even 
smaller radius, whereas AV18 plus Illinois-2 leads to 
a radius [H) that is closer to the experimental value. 
Note however that the gs quadrupole moment obtained 
with that interaction is significantly larger than experi- 
ment [3l| . 

Also the calculated magnetic moment and reduced Ml 
transition probabilities with JISP16 compare well with 
available experimental data for all the cases shown in 
Table H 

The (0 + , 1) excited state of 6 Li is the isobaric analog of 
the ground state of 6 He, and can be used to calculate the 
Gamow-Teller transition between 6 Li and 6 Hc. Assum- 
ing isospin symmetry, the Gamow-Teller matrix element 
Mqt is to good approximation given by 



A I, 



GT 



(21) 



and is related to the half-life [26J through 



|M GT | 2 



1 2^ 3 ln2/(G 2 |K d 



fv 9A 



(fT 1/2 )tml 



(22) 



1.2695(29) is the axial constant, — 



fv 



where gA 

1.00529 accounts for the difference in the statistical rate 
function of the vector and axial- vector transitions, m e is 
the mass of the electron, G = 1.166371(6) • 10- n MeV" 2 
is the Fermi coupling constant, and V u d = 0.9738(4) 
is the CKM matrix element that mixes the quarks in- 
volved in the decay Our Mqt result, presented in Ta- 
ble HI compares quite well to that calculated in [26| using 
the hyperspherical-harmonic expansion method with the 
same (JISP16) interaction; they also obtained a value 
of M GT = 2.227. It is interesting to note that Ref. HI 
found the exchange current corrections to the GT matrix 
element to be of the order of a few percent. Ref. [26| also 
presents the extrapolated gs binding energy of 6 Li with 
JISP16 as 31.46(5) McV which is in excellent agreement 
of our 31.49(3) McV result presented in Table [D 



B. Density distributions 

A closer look at the 3-dimensional one-body densities, 
free of spurious cm effects, would help to develop a phys- 
ical intuition for the ab initio structure of a nucleus. Al- 
though it is easier to perform the deconvolution of the 




FIG. 5: (Color online) The y = slice of the translationally- 
invariant matter density in the x-z plane for the gs of 6 Li (top, 
J = 1) is contrasted with the density for the first excited 
state (bottom, J — 3). These densities were calculated at 
jV max = 16 and TiQ, = 12.5 MeV. 



cm density after integrating out all angle-dependence one 
can also deconvolute the full 3-dimcnsional density (see 
the appendix for details). However, a detailed investi- 
gation of the numerical convergence is impractical for 
these 3-dimensional densities. We therefore present all 
our 3-dimensional density distributions in the largest ba- 
sis space at Ml = 12.5 MeV only, where the RMS radius, 
as well as the quadrupole moments (which are closely 
related to the shape of the wavefunction) , appear to be 
reasonably converged. 

In order to produce the density that represents the 
actual shape of a specific state of a nucleus in a 
translationally-invariant (inertial) frame, we set Mj = J 
for all our calculations of the local density. That is, we se- 
lect the maximal positive angular momentum projection 
along the axis of quantization, the z-axis. This seems 
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like the natural choice since the quadrupole moment Q is 
defined as the E2 matrix element at J = Mj (or equiv- 
alcntly, to the reduced E2 matrix element). Note that 
even though we calculate 3-dimensional density distribu- 
tions, our results are symmetric under rotations around 
the z-axis: the wavefunctions have azimuthal symmetry. 

Fig. [5] shows the matter density for the lowest two 
states of 6 Li. Both states are oblate, though the (l + ,0) 
state is nearly spherical whereas the (3 + ,0) state is 
strongly oblate. Indeed, the relative deformation of the 
translationally-invariant densities for the gs and first ex- 
cited state is implied by the results in Table Q] for their 
respective quadrupole moments. The gs has a negative 
calculated quadrupole that is near zero, in close agree- 
ment with experiment. On the other hand, the first ex- 
cited state has a large negative calculated quadrupole 
moment. 

It is worth commenting that our use of the terms "pro- 
late" and "oblate" characterize the shapes in the inertial 
frame, not a body-fixed axis as is common for discus- 
sions of shapes in the collective model In the in- 
ertial frame of reference positive quadrupole moments 
correspond to prolate shapes and negative quadrupole 
moments correspond to oblate shapes. 

Fig. |5] illustrates the effect of Mj on the density dis- 
tribution as we see the oblate shape of the density at 
Mj = 3 (top panel) morph into the prolate shape at 
Mj = (bottom panel). Calculating the density when 
J =/= AIj gives a density whose azimuthal symmetry axis 
is not aligned with the spin. E.g. at Mj = the symme- 
try axis lies in the x-y plane, perpendicular to the z-axis. 
The oblate shape we found at Mj = 3 is now also per- 
pendicular to the z-axis. On the other hand, we do have 
an azimuthal symmetry around the z-axis. Therefore, 
what we obtain is an oblate shape perpendicular to the 
z-axis, but with its principle axis (symmetry axis) av- 
eraged over all directions in the x-y plane. This results 
in a prolate shape at Mj = 0, as we see in the bottom 
panel of Fig. [5] Note that we see the same for e.g. the E2 
matrix element: with Mj = 0, the E2 matrix element for 
this state is positive, but the corresponding quadrupole 
moment, shown in Table HI is negative (and independent 

Of Mj). 

C. 7 Li 

For 7 Li we evaluated the low-lying states in basis 
spaces up to N max = 14. We only consider isospin | 
states - the lowest isospin | has more than 10 MeV 
excitation energy. The lowest five states in the exci- 
tation energy spectrum compare well with experiment 
and the correct level ordering is preserved as shown in 
Fig. U\ The excitation energy of four of these five states 
shows rapid numerical convergence with iV max and stabil- 
ity with respect to variations in the HO energy. However, 
the convergence of the lower of the two | states is sig- 
nificantly slower. Indeed, experimentally this state has a 




x (fm) 



FIG. 6: (Color online) The y = slice of the translationally- 
invariant matter density in the x-z plane for first excited 3 + 
state of 6 Li with Mj = 3, 2, 1, from top to bottom. 



10 



10 



8 



0) 



M 

Si 

0) 
3 



'Li 

5/2' Exp. 
5/2- Exp. 



7/2" Exp. 



1/2" Exp. 



N =8 

" max 

N =10 

max 

N =12 

max 

N =14 

" max 



12.5 



15 17.5 20 

fiQ, (MeV) 



22.5 25 27.5 



FIG. 7: (Color online) Excitation energies for selected excited 
states of 7 Li are shown as a function of Ml at various V max 
values as indicated in the legend, and compared to experi- 
ment. 



large width of 0.88 MeV, whereas the width of the other 
states is less than 0.1 MeV. Thus, as in 6 Li, we again 
observe a good correlation between experimental width 
and convergence rate of excitation energies. 

Our results with JISP16 for selected spectral and other 
observables of 7 Li are summarized in Table |TT] and com- 
pared with experiment when available. The error esti- 
mates for the excitation energies are calculated as dis- 
cussed above. From this table we see that the gs energy 
is underbound by about 0.67 MeV. The gs energy and 
excitation energies compare favorably to other methods 
and interactions [3(1 H3j- For the RMS radius and the 
quadrupole observables (which are not converged) we list 
our results from the largest basis space, iV max = 14, at 
Ml = 12.5 MeV; they are in reasonable agreement with 
the available experimental data, given the basis space 
limitations and long-range nature of these operators. 

Fig. [8] displays the magnetic dipole moments and re- 
duced magnetic dipole transition probabilities to the gs 
for selected states of 7 Li. These observables converge 
quickly, and for calculations at -/V max = 14, are almost 
independent the HO energy In fact, most of the mag- 
netic observables are already reasonably well converged 
(to within 10%) at N max = 8, with the noticable excep- 
tion of the two | excited states. This is partially due to 
a strong state mixing between these two states. We re- 
quire larger basis spaces to fully differentiate these states, 
because they are close together in energy, and their quan- 
tum numbers are identical. 

Our estimate for the infinite basis space results for 
magnetic dipole observables is based on the residual de- 
pendence on N max and TiVL over a 10 MeV window in HQ. 
This window does include the optimal h£l for the ex- 
trapolations and the variational upperbound, but is not 
necessarily centered around these values. Our numerical 
error estimate is the RMS sum of the variation with Ml 



7 Li 


Expt. 


JISP16 


AV18/IL2 


CD-B 


INOY 




39.244 


38.57(4) 


38.9(1) 


35.56 


39.6(4) 


(r 2 ) 1/2 

V ppi 


2.30(5) 


2.2 


2.25(1) 


2.22 


2.05 


Mr) 


0.477 


0.52(6) 


0.2(1) 








4.630(1) 


5.25(5) 


4.9(1) 








6.680(50) 


7.1(2) 


6.6(1) 








7.460(10) 


8.1(1) 


7.2(1) 






Q(§~) 


-4.06(8) 


-3.2 


-3.6(1) 


-3.20 


-2.79 




- 


-5.0 








cdD 


- 


-6.0 










- 


2.3 










3.256 


2.954(5) 


3.168(13) 


3.01 


3.02 




- 


-0.76(1) 










- 


3.3(1) 








Mir) 


- 


-0.90(2) 










- 


-0.39(5) 








B(E2;|-) 


15.7(10) 


10.2 


16.2(5) 






B(E2;|") 


3.4 


5.1 


9.92(14) 






B(E2;| 1 ) 


- 


1.5 








B(E2;| 2 ) 




< 0.1 








B(M1;§1 


4.92(25) 


3.89(2) 


4.92(7) 






B(M1;§-) 




0.002(1) 








B(M1;§-) 




0.02(1) 









TABLE II: Selected 7 Li observables calculated up through 
-Vmax = 14, with the same units as in Tabled] The energies are 
obtained from extrapolations to the infinite basis space, and 
the magnetic dipole observables are nearly converged, with 
error estimates as discussed in the text; the RMS point-proton 
radius and electric quadrupole observables are evaluated at 
Ml = 12.5 MeV. Experimental values are from Refs. [2l|, [2?], 
AV18/IL2 results are from Refs. [Tl],[lHil and include 
meson-exchange corrections for the dipole observables; CD- 
Bonn and INOY results are from Ref. [32] , and were calculated 
at A r maa; =12 and ^,£1=11 and 16 MeV respectively for CD- 
Bonn and INOY, with the INOY gs energy extrapolated to 
the infinite basis space. 



over this window and the difference between the results 
in the two largests N max calculations (rounded up), i.e. 
treating the variation with each of the two basis space 
parameters as independent sources of numerical uncer- 
tainties. 

With JISP16, the magnetic moment of the gs is about 
10% too low compared to experiment. This could eas- 
ily arise from our neglect of meson-exchange currents in 
our current calculations. GFMC calculations with AV18 
plus Illinois-2 three-body forces [ll[ found that the mag- 
netic moment of the 7 Li gs receives a 10% correction 
from meson-exchange currents, changing the magnetic 
moment from 2.9 [In to 3.2 //jy. It is quite remarkable 
that our result for the magnetic moment, with the naive 
pointlikc Ml operator, is in fact quite close to the results 
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FIG. 8: (Color online) Dipole moment (top) and B(M1) tran- 
sition to the gs (bottom) as a function of the basis htt for 
selected excited states of 7 Li at various A^ max values as indi- 
cated in the legend. The experimental value for the B(M1) 
transition from the excited ^ to the ground state is indi- 
cated on the left side of the lower panel and is about 20% 
larger than the theoretical result. 




FIG. 9: (Color online) The y = slice of the gs matter den- 
sity of 7 Li before (top) and after (bottom) deconvolution of 
the spurious cm motion. These densities were calculated at 
jV max = 14 and HQ. = 12.5 MeV. 



obtained with CD-Bonn, with INOY, and with AV18 plus 
Illinois-2, all about 10% below the experimental datum. 
Apparently, this observable is not very sensitive to the 
details of the interaction. (Note that the exchange cur- 
rent correction to the 6 Li gs magnetic moment was only 
2% in Ref. [ll|.) 

The B(M1) from the \~ to the gs is about 20% too 
low compared to experiment. Again, this is in qualitative 
agreement with the findings of Ref. [ll|: with AV18 plus 
Illinois-2 there is about 10% increase in the Ml transi- 
tion matrix element due meson-exchange currents, which 
results in a 20% increase in the corresponding B(M1). 

The effect of the cm motion on the density is shown 
in Fig. [JJ] for the gs of 7 Li. The top panel shows the 
space-fixed (sf) density including the cm motion, Pgf(r), 
whereas the bottom panel shows the translationally- 



invariant density, Pt\(f). The smearing of the den- 
sity due to the cm motion leads to a diminished cen- 
tral density; the sf density has a central value of 0.204 
nucleons/fm 3 while the ti density has a central value of 
0.233 nucleons/fm 3 . 

The cm motion smearing spreads out the sf density 
leading to a slower falloff and a larger radius than the 
ti density. Furthermore, the ti density has a more pro- 
nounced oblate shape than the sf density, as would be ex- 
pected from smearing with a spherically-symmetric func- 
tion that averages out the non-spherical details. In order 
to characterize the degree of deformation, we compare 
the ratio of the long axis to the short axis of the ellip- 
tical density slices. The ratio of the long to short axes 
at half central density is 1.78 for the ti density and 1.60 
for the sf density. Note that the extent of the smearing 
effect from cm motion depends on the HO energy of the 
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FIG. 10: (Color online) The y = slices of the translationally- 
invariant proton densities for the first excited | state (top) 

and for the second excited |~ state (bottom) of 7 Li. These 
densities were calculated at -/V max = 14 and HQ, — 12.5 MeV. 



basis. The sf density depends on TiQ,, even in the limit 
iV max — !> oo, whereas the ti density becomes independent 
of the basis in this limit. 

Fig.[TOJcontrasts proton densities of the fourth and fifth 
excited states of 7 Li. Although their quantum numbers 
are equal, (J*,T) = (| , h), they have other features 
that make them quite distinct. Experimentally, the first 
excited § is broad, whereas the second excited | is 
narrow: their widths are 0.88 MeV and 0.09 MeV re- 
spectively. Indeed, our calculated excitation energy is 
better converged for the higher of these two states. Fur- 
thermore, our calculations show significant differences in 
their structure: the first excited § has a large negative 
quadrupolc moment, whereas the second has a moder- 
ate positive quadrupolc moment (see Table ITT|) . Indeed, 
the density in the top panel of Fig. [TO] is strongly oblate, 



8 Li 


Expt. 


JISP16 


AV18/IL2 


CD-B 


E b (2+) 


41.277 


40.3(2) 


41.9(2) 


35.82 


( r pp) 


2.21(6) 


2.1 


2.09(1) 


2.17 


£x(l + ) 


0.981 


1.5(2) 


1.4(3) 




£ x (3+) 


2.255(3) 


2.8(1) 


2.5(3) 




£x(4+) 


6.53(2) 


7.0(3) 


7.2(3) 




Q(2 + ) 


3.27(6) 


2.6 


3.2(1) 


2.78 


Q(i + ) 


- 


1.2 






Q(3 + ) 


- 


-2.0 






Q(4 + ) 


- 


-3.4 






M(2+) 


1.654 


1.3(1) 


1.65(1) 


1.24 


m(i + ) 


- 


-2.2(2) 






A*(3+) 


- 


2.0(1) 






A*(4+) 


- 


1.84(1) 






B(E2;1+) 


- 


1.9 






B(E2;3+) 




4.6 






B(E2;4+) 




1.9 






B(M1;1+) 


5.0(16) 


3.7(2) 






B(M1;3+) 


0.52(23) 


0.25(5) 







TABLE III: Selected 8 Li observables calculated up through 
-Vmax = 12, with the same units as in Table [I] The ener- 
gies are obtained from extrapolations to the infinite basis 
space, and the magnetic dipole observables are nearly con- 
verged, with error estimates as discussed in the text; the RMS 
point-proton radius and electric quadrupole observables are 
evaluated at hi} — 12.5 MeV. Experimental values are from 
Refs. [H], [H, Hi|- AV18/IL2 results are from Refs. [U, EH 
and does not include meson- exchange corrections for the mag- 
netic moment; CD-Bonn and INOY resuts are from Ref. [32 |. 
and were calculated at N max — 12 and HQ— 12 and 16 MeV re- 
spectively for CD-Bonn and INOY, with the INOY gs energy 
extrapolated to the infinite basis space. 

whereas the bottom panel shows a moderately prolate 
shape (we note again that the densities are symmetric 
around the azimuthal axis, which is the vertical axis is 
these plots). Another noteworthy difference is observed 
in the magnitude of the central proton density: the more 
diffuse (| , i)i state has a central proton density of only 

0.08 protons/fm 3 while the (| , state has a central 
proton density of 0.12 protons/fm 3 , 50% higher. 

D. 8 Li 

Table IIIII presents a capsule view of selected spectral 
and other observables for 8 Li calculated with JISP16 
and compared with experiment when available. The er- 
ror estimates for the excitation energies and magnetic 
dipole observables are calculated as discussed above. For 
the RMS point-proton radius and the charge quadrupolc 
observables (which are not converged) we list our re- 
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FIG. 11: (Color online) Excitation energies for select excited 
states of 8 Li are shown as a function of Ml at various iV ma x 
values as indicated in the legend. 



suits from the largest basis space, -/V max = 12, at hfl = 
12.5 MeV. From this table we see that the gs energy is 
underbound by about 1.0 MeV. 

In addition to the gs, we examined several narrow low- 
lying states: the lowest two excited states, which are 
narrow states, with a width of 33 keV or less, as well as 
a narrow low- lying 4 + state at 6.53 MeV with a width 
of 35 keV. We do not consider isospin 2 states - the low- 
est isospin 2 has more than 10 MeV excitation energy. 
The excitation energies obtained with JISP16 compare 
reasonably well with the experimental excitation ener- 
gies, though the level splittings are a bit too large (see 
Fig. [TTj) . The convergence of the spectrum is similar as 
for the other Li-isotopes, though the convergence of the 
4+ state is somewhat slower than expected based on its 
small width. In our calculations there are several addi- 
tional states below this 4 + state, which is typically the 
8th state in our calculated spectrum, depending on the 
basis space; experimentally, it is the 7th observed state. 

The convergence properties of the calculated magnetic 
dipolc obscrvables are similar to those for 6 Li and 7 Li dis- 
cussed above. Because we can only go up to iV max = 12, 
the numerical error estimates are slightly larger than for 
6 Li and Li. The gs magnetic moment is approximately 
20% lower than experiment. This seems reasonable, in 
light of our discussions of magnetic moments above, that 
this discrepancy is at least partially due to the fact that 
we do not incorporate meson-exchange currents. 

Note that CD-Bonn gives a similar magnetic moment 
as JISP16, but that INOY provides a moment that is 
slightly closer to experiment, and AV18/IL2 gives a mag- 
netic moment in excellent agreement with data; however, 
the contributions of meson-exchange currents were not 
evaluated for this magnetic moment. 

The reduced B(M1) transition probabilities from the 
1 + and 3 + to the gs are 20% and 50% lower than ex- 



periment, but the experimental error bars are large. The 
magnetic moment of the 4 + state is remarkably well con- 
verged, despite its excitation energy not being very well 
converged. 

The quadrupolc moments and reduced B(E2) transi- 
tion probabilities are not well converged due to basis 
space limitations, as discussed above. In spite of these 
limitations, the electric quadrupole moments allow us to 
qualitatively understand the shape of the proton densi- 
ties of these states: prolate for the 2 + gs and the first 
excited 1 + , but oblate for the first excited 3 + and 4 + 
states. Indeed, that is what we see in the left-hand pan- 
els of Fig. [T^l where we plot these densities for the gs and 
for the 4+ state. 

Interestingly, the neutron density differs by more than 
a simple scale change from the proton density for these 
two states, as can be seen from the right-hand panels of 
Fig. [T2] In the 2 + state, the deformation of the neutrons 
is significantly larger than that of the protons, whereas 
in the 4 + state, the deformation of the neutrons is much 
smaller than that of the protons. 

A case of special interest can be seen in the top right 
panel of Fig. [T^l or in more detail in Fig. Q2J In the left 
panel of Fig.[L|]wc clearly see non-trivial neutron cluster- 
ing that is obfuscated in the sf frame (right), highlighting 
the importance of the deconvolution procedure and the 
significance of the ti density. Furthermore, the ti density 
(left panel) has a significantly higher density in the cen- 
tral region than the sf density (right panel). Both the ti 
and the sf densities are normalized to give the same in- 
tegrated density of five neutrons, so that means that the 
sf density is smeared out over a larger region, and falls 
off to zero slower than the ti density. This is also evident 
in Fig. O where we contrasted the sf and ti matter den- 
sity of the gs of 7 Li: the ti central density is significantly 
higher than the sf central density. 

Another way of visualizing these ti densities is by plot- 
ting their multipolc components, p( K ^(r) (we omit the 
"ti" subscript for compactness of notation). In order to 
exhibit their long distance features, we present the mag- 
nitudes of the p( K \r) in Fig. [T3]for the 2 + gs and the 
first excited 4 + state of 8 Li as semilog plots out to a ra- 
dius of 6 fm. To illustrate the radial and angular nodal 
structure together out to 4 fm, we present in Fig. [T5] 
the y = slices of the multipolc components for the gs 
proton density (left panels) and gs neutron density (right 
panels) defined as their full contribution to the total den- 
sity. That is, Fig. [T5l displays slices of the respective full 
terms contributing to the sum given in Eq. (|10[) . Hence, 
the sum of the proton (neutron) multipoles in Fig. 1151 
produces the top left (top right) panel of Fig. out to 
4 fm. 

Let us consider Fig. [TJ] in more detail. Qualitatively, 
the multipole components look very similar for the gs pro- 
tons and gs neutrons. The main difference seems to be 
that the proton densities fall off more rapidly with r than 
the neutron densities for all three multipole components. 
This is understandable since this is a neutron-rich system 
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FIG. 12: (Color online) The y — slice of the translationally-invariant proton (left) and neutron (right) densities of the 2 + gs 
(top) and the first excited 4 + state (bottom) of 8 Li. These densities were calculated at -/V max = 12 and htt = 12.5 MeV. 



and the single-neutron removal energy is less than that 
of the single-proton removal energy. Note however that 
the clustering of the neutrons in the gs of 8 Li (left panel 
of Fig. I13p is not evident from the multipole components 



P 



of the neutron density displayed in the top right 
panel of Fig. Q31 That is, even though the radial multi- 
pole densities of the protons and neutrons look qualita- 
tively similar, the corresponding 3-dimcnsional densities 
look qualitatively different. Apart from overall differ- 
ence in the asymptotic behavior, the biggest difference 
between the proton and neutron multipole densities is in 
the interior of the nucleus, below 2.5 fm, where the mag- 
nitude of the proton hexadecapole density is up to 60% 
larger than that of the neutron hexadecapole density. On 
the other hand, in the exterior region, beyond 4 fm, the 
neutron hexadecapole density is more than an order of 
magnitude larger than the proton density. 

The monopolc proton and neutron densities of the first 
excited 4 + state of 8 Li seen in Fig. [THare similar to those 



of the ground state, with the proton density falling off 
more rapidly than the neutron density. On the other 
hand, the higher multipole components p( K ^(r) of the 
first excited 4 + state of 8 Li look qualitatively quite dif- 
ferent than those of the ground state. In addition, the 
protons and the neutrons of the excited 4 + state differ 
from each other in their higher multipole components. 
The quadrupolc density of the neutrons in the 4 + state 
has a node, in contrast to that of the protons. Neither the 
proton nor the neutron hexadecapole density has a node 
in the 4 + state, whereas both the proton and the neutron 
hexadecapole densities have a node in the ground state. 

We turn our attention now to the nodal structures of 
the multipole components of the charge-dependent den- 
sity distributions for the gs of 8 Li in Fig. [15] Since each 
panel represents a term contributing to Eq.( I10p. we ob- 
serve the angular nodal structure governed by the respec- 
tive spherical harmonic factor. The radial nodes arise 
from the radial functions p( K >(r) displayed in Fig. [TJ] 
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x (fm) x (fm) 

FIG. 13: (Color online) The y = slice of the translationally-invariant neutron density (left) of the 2 + gs of 8 Li. The space-fixed 
neutron density for the same state is on the right. The contour labels give the density in neutrons/fm 3 . These densities were 
calculated at A max = 12 and HQ = 12.5 MeV. 



The significant differences between the proton and the 
neutron magnitudes at each multipolarity reflect the dif- 
ferences in the magnitudes of the radial distributions of 
the multipoles seen above in the top two panels of Fig.lT4"l 



From both Figs. [TH and Q~5l one observes that the gs 
proton monopole density has a slightly higher central 
value than the neutron gs monopole density, but it falls 
off much more rapidly with r than the neutron density. 
On the other hand, the gs neutron quadrupolc density 
has larger (in magnitude) features than the gs proton 
density, both at small and large distances. And again, 
the gs hexadecapole densities have the most interesting 
features, including the nodal structure at about 3 fm. In 
the interior region of the nucleus, the gs proton hexade- 
capole density has more pronounced features, whereas in 
region beyond 3 fm the gs neutron hexadecapole den- 
sity has more pronounced features. However, none of 
these multipole density plots show any hint of clustering 
of the neutrons; nevertheless, the plots of the three gs 
neutron multipole densities of Fig. [15] add up to give the 
total translationally-invariant density of the left panel of 
Fig. [UH which does indicate two neutron clusters, with 
their centers separated by about 1.5 fm. 



Although these densities are not (yet) fully converged, 
we feel that the qualitative features will persist in the 
limit of a complete basis. In particular salient differences 
between different states and/or between the proton and 
neutron densities arc likely to survive in that limit. 



IV. SUMMARY AND OUTLOOK 

We have performed no-core full configuration calcula- 
tions for the Lithium isotopes, 6 Li, 7 Li, and 8 Li with the 
realistic NN interaction JISP16. Several observables ob- 
tained (gs energies, excitation energies, magnetic dipole 
moments and reduced magnetic dipole transition prob- 
abilities) compare well with both experiment and alter- 
nate methods and interactions. For certain observables 
that are more sensitive to long-range correlations (the 
RMS radius, electric quadrupolc moments, and reduced 
quadrupole transition probabilities) we were unable to 
obtain full convergence, though they also compare favor- 
ably with alternate methods and interactions. 

One-dimensional and three-dimensional 

translationally-invariant one-body density distribu- 
tions were calculated for various ground and excited 
states of 6 Li, 7 Li, and 8 Li. These one-body density 
distributions provide an excellent framework for visu- 
alization of nuclear shape distortions and clustering 
effects. The associated one-body density matrix in the 
HO basis provides a compact form of all of the quantum 
one-body information for a given nuclear state. 

To improve our convergence, especially for matrix el- 
ements of long-range operators, we would require sig- 
nificant increases in basis space sizes (increased N max ) 
and/or alternatives to the HO single-particle basis. Re- 
cent advances in the "Importance- Truncated No-Core 
Shell Model" ESIU, the "Symmetry- Adapted No-Core 
Shell Moder_p| and the "No-Core Monte Carlo Shell 
Model" [43|, [44 1 are promising new methods for accessing 
much larger basis spaces. 
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FIG. 14: (Color online) The multipole components p ti (r) of the proton (left) and neutron (right) densities of the 2 + gs (top) 
and the first excited 4 + state (bottom) of s Li. These densities were calculated at iV max = 12 and HQ, = 12.5 MeV. Monopole 
and quadrupole distributions for the gs are all positive. The K = 4 distributions for the gs are negative in the interior and 
positive in the tail region. For the 4 + state, the monopoles are positive while the quadrupole is negative for the protons and 
negative for the interior of the neutrons. Both K = 4 distributions are positive for the 4 + state. 



Further advances in NN interactions, as well as three- 
body forces, could also help resolve some of the residual 
differences between theory and experiment. Of course, 
there is also the possibility that four-body forces may 
play a significant role. 
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Appendix A: Deconvolution of the one-body density 

In our many-body framework, we use the Slater de- 
terminants of a single-particle Harmonic Oscillator (HO) 
basis. In this basis, we can define a space-fixed One-Body 
Density Matrix (OBDM) by its matrix elements 

pg = (* f \atap\%), (Al) 

where a and /3 stand for a set of single-particle quantum 
numbers (n a , l a , j a , m a , t Zi01 ) and {np,lp,jp,mp,T Z! p). 
These matrix elements p^ a , together with the expressions 
for the single-particle wavefunctions ip a (f), completely 
determine the OBDM in coordinate space 

pV.O = Y,P^^(r)Mn- (A2) 

a,0 
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x (fm) x (fm) 



FIG. 15: (Color online) The y = slices of the translationally-invariant proton (left) and neutron (right) densities of the 2 + 
gs of 8 Li. From top to bottom, we present the monopole, quadrupole and hexadecapole densities respectively. These densities 
were calculated at -/V max = 12 and %Q, — 12.5 MeV. 
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The local one-body density becomes 



P*t 



(A3) 



where 



< f\a > 

E (l a mis a m s \j a i 



"j<x/ TTia l a rni 



(r) 



mi ,171 s 



(A4) 



with 4>n a i a mi{'f f ) as the usual three-dimensional HO 
coordinate-space wavefunction and Xs a m a is the Pauli 
spinor for the nuclcon. While we retain the same basis 
functions ip a for neutrons and for protons (subject to po- 
tentially different single-particle state cutoffs as dictated 
by -/V max ), our many-body states treat the neutrons and 
protons independently so isospin is not a conserved quan- 
tity. The one-body density distribution is normalized to 
the number of nucleons 



J &{r)d\ = A. 



(A5) 



This local one-body density distribution includes con- 
tributions from the cm (center of mass) motion of the 
many-body wavefunctions ^(fi Va), hence the sub- 
script sf (space- fixed) . However, because of the ex- 
act factorization of the cm wavefunction and the ti 
(translationally-invariant) wavefunction, the sf density 
can be expressed as a convolution of the ti density dis- 
tribution p t i with the cm density distribution p cm via 



PS(0 = 



Pti (f-R) p? m (R) d 3 R 



(A6) 



For a HO basis, is a simple Gaussian (the gs den- 
sity of -Hem) that smears out any non-spherical features 
of pti- This smearing also introduces into p an unde- 
sired dependence on the basis parameter TiVl that could 
mask the convergence of pti- Given our expansion in a 
HO basis, we have analytically evaluated Eq. (A7) using 
our OBDMs with Mathematica (45| . This produces an 
analytic expression for the ti OBDM. 

In order to display properties of the density pti, as well 
as investigate its convergence, we deconvolute the spuri- 
ous cm density using well-known Fourier methods [l5j |. 
The ti density is given by 



Pti(r) 



F~ 



F[p<? m (R)} 



(A7) 



where F[f(r)] is the 3-dimensional Fourier transform of 

For an ^4-nucleon eigenstate with total angular mo- 
mentum J and projection M, and possibly additional 
quantum numbers A, denoted by \AXJM), we can evalu- 
ate the local space-fixed density by evaluating the matrix 



element [46| (henceforward, we omit the superscript "ST 
for compactness of notation) 



p sf (f) = (AXJM\p si (r)\AXJM) 
using the local density operator 
A 



(A8) 



Psf(r) = ^<5 3 (r-r fc ) 



k=l 
A 



E — 2^ Y,Yi* m (?k)Yr(r) (A9) 



k=l 



where f is the unit vector in the direction r, and YJ m (f) 
is a spherical harmonic. Note that it has the property 



Y- m (r) = (-l)T m (r). 



(A10) 



To efficiently perform the deconvolution of the spurious 
cm density and the ti density, we make a multipole ex- 
pansion of the local density [46[ 



Psi(r) = E 



(JMK0\JM) {K) 



K 



V2J+1 



Y£>(r)pl?'(r) (All) 



where p^ (r) is the K th multipole of the sf density. The 
same equation with "sf ' replaced by "ti" provides the cor- 
responding expansion for the local ti density. For initial 
and final states with spin J, the multipoles range from 
K = to K = 2J. This multipole expansion greatly 
simplifies the Fourier transforms needed for the deconvo- 
lution, as we will see at the end of this appendix. 

With a HO single-particle basis, each multipole is given 
by 



Psf 



(r) = ^RmiA^Rn^ir) -^{h\ji\\Y K \\l 2 \j2) 
x{AXJ\\(al ihji ~a n2l2j2 nAXJ), (A12) 



where K = \J2K + 1. The i?„i(r)'s are the radial comp- 
nents of the HO wavefunction 



Rni(r) 



2(2i/) z+3 / 2 r(n + 1) 



T(n 



1/2 



r2 L l f*{2vr 2 ) 



(A13) 



with L n 2 the associated Laguerre polynomials and v = 
mc 2 Ml/ (2h 2 c 2 ). The reduced matrix element of a spher- 
ical harmonic in Eq. (|A12[) can be written as 

1 



(h^ji\\Y K \\l 2 ^j2) = 



x(-iy* + i (h0h0\K0) I f 2 j * f I (A14) 

using a Wigner-6J coefficient. (We use the Conden- 
Shortley convention for the Clebsch-Gordan coefficients, 
defined in Ref. 0.) 
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Finally, (AXJ\\(al ihn ~a n2l232 )^\\AXJ) in Eq. flXT5J 
represents a reduced matrix element of the K th multipolc 
of the OBDM. The K th multipolc operator for initial and 
final states with the same Mj can be written as 

(ai lhH an 2 i 2n ) (K) = ^(-l)^ +m Mii^i2-m,|irO) 



X a n 1 hj 1 m j an 2hj2m j 



(A15) 



Note that in Eq. (|Al2j) we use Mj -independent reduced 
matrix elements, defined by the Wigncr-Eckart theorem. 
For a generic operator Tk k the reduced matrix element 
is defined by 



(XfJfWTxWXiJj = J f 



(\fJfMf\TKk\XiJiMi\ 
(J f M f Kk\JiMi) 



(A16) 



provided that the Clcbsch-Gordan coefficient in the de- 
nominator is not zero. Thus the reduced matrix elements 
of the K th multipole of the OBDM are given by 

(AXJ\\(ai ihn ~a n2l2j2 )^\\AXJ) = 

= ^ 2J + l Y(-lY 2+m i (jim 7 7 2 - mAKO) 

rrij 

x {AXJM\a{ ihjim] a n2l2nm] \AXJM) . (A17) 

where (AXJM\al ilijim a n2 i 2 j 2mj \AXJM) arc the OB- 
DME's in HO spa ce defined in Eq. (TAT]) . Note that 
the factor \J2J + l/( JMK0\ JM) in this expression can- 
cels against the same factor in the expression for p s f(f), 
Eq. (TATTj) . 

We can now efficiently perform the deconvolution, by 
using the Fourier transform properties of the multipole 
expansion of the local density 



(A18) 



where the multipole component of the density in momen- 
tum space is expressed as 

pifHq) = 4tt / j K (qr) p(f>(r) r 2 dr (A19) 



with jx the spherical Bcssel Functions of the first kind. 
Thus the deconvolution of each multipole gives 



where 



872 3/2 f°° e-^ s 



qR 



R dR 



(A21) 

For spherically symmetric nuclei, this deconvolution 
simplifies even further because we only have one term in 
the multipole expansion, K = 



P& 0) (r) 



sin(qr) pf { } (q) 
2tt 2 J qr p C m{q) 



q 2 dq, (A22) 



and the 3-dimcnsional ti density is simply 



Pti(f) 



2V^F 



(A23) 



without any angular dependence. 

Another advantage of the multipole expansion is that 
it allows for a straightforward calculation of the (sf or 
ti) density for any magnetic projection M, once the mul- 
tipolcs p( K >(r) are known. The multipoles p^ K \r) arc 
completely determined from reduced matrix elements, 
wich do not depend M. The only dependence of p s f{r) 
on M is entirely through the explicitly M-dependent 
Clebsch-Gordan coefficients in Eq. (|A11I) . 

Finally, the expressions here are specific for the local 
static density (the same final and initial state), but the 
extension to local transition densities is straightforward. 



[1] A. M. Shirokov, J. P. Vary, A. I. Mazur and T. A. Weber, 

Phys. Lett. B 644, 33 (2007) 
[2] A. M. Shirokov, J. P. Vary, A. I. Mazur, S. A. Zaytsev 

and T. A. Weber, Phys. Lett. B 621, 96 (2005) 
[3] A. M. Shirokov, A. I. Mazur, S. A. Zaytsev, J. P. Vary 

and T. A. Weber, Phys. Rev. C 70, 044005 (2004) 
[4] P. Maris, A. M. Shirokov and J. P. Vary, Phys. Rev. C 

81, 021301 (2010) 



[5] P. Maris, J. P. Vary and A. M. Shirokov, Phys. Rev. C 

79, 014308 (2009) 
[6] S. K. Bogner, R. J. Furnstahl, P. Maris, R. J. Perry, 

A. Schwenk, J. P. Vary, Nucl. Phys. A 801, 21-42 (2008) 
[7] B. D. Keister, W. N. Polyzou, J. Comput. Phys. 134, 

231-235 (1997) 

[8] M. A. C aprio, P. Maris an d J. P. Vary, Phys. Rev. C. (to 
appear) ; I arXiv: 1208.41561 



20 



P. Navratil, J. P. Vary and B. R. Barrett, Phys. Rev. C 
62, 054311 (2000). 
K. Suzuki and S. Y. Lee, Prog. Theor. Phys. 64, 2091 [29 
(1980); K. Suzuki, ibid. 68, 246 (1982); 68, 1999 (1982); 
K. Suzuki and R. Okamoto, Prog. Theor. Phys. 92, 1045 [30 
(1994). 

L. E. Marcucci, M. Pervin, S. C. Pieper, R. Schiavilla, [31 

R. B. Wiringa, Phys. Rev. C 78, 065501, (2008). 

D.H. Gloekner and D.R. Lawson, Phys. Lett. B 53, 313 [32 

(1974). 

P. Sternberg, E. G. Ng, C. Yang, P. Maris, J. P. Vary, [33 
M. Sosonkina, and H. V. Le, in the Proceedings of the 
2008 ACM/IEEE conference on Supercomputing IEEE [34 
Press, Piscataway, NJ, 15:1-15:12 (2008). [35 
P. Maris, M. Sosonkina, J. P. Vary, E. G. Ng, and 

C. Yang, Procedia CS 1, 97-106 (2010). 
P. Morse and H. Feshbach, Methods of Theoretical [36 
Physics, McGraw-Hill, New York (1953). 

B. G. Giraud, Phys. Rev. C 77, 014311 (2008). [37 
P. Navratil, Phys. Rev. C 70, 014317 (2004). 
I. Talmi, Simple Models of Complex Nuclei, Harwood [38 
Academic Publishers, Switzerland (1993). 
A. Bohr and B. Mottelson, Nuclear Structure, W.A. Ben- 
jamin, New York (1969). [39 
J. L. Friar, J. Martorell, and D. W. L. Sprung, Phys. 
Rev. A 56, 4579 (1997). [40 
W. Nortershauser, T. Neff, R. Sanchez, I. Sick, Phys. [41 
Rev. C 84, 024307 (2011). [42 
R. Sanchez et al, Phys. Rev. Lett. 96, 033002 (2006). 
A. M. Shirokov, A. I. Mazur, J. P. Vary and E. A. Mazur, [43 
Phys. Rev. C 79, 014610 (2009) 

I. Stetcu, B. R. Barrett, P. Navratil and J. P. Vary, Phys. 
Rev. C 73, 037307 (2006) [44 
R. Roth, T. Neff, H. Feldmeier, Prog. Part. Nucl. Phys. 
65, 50-93 (2010) [45 
S. Vaintraub, N. Barnea, and D. Gazit, Phys. Rev. C 79, 
065501 (2009) [46 
F. Ajzenberg-Selove, Nucl. Phys. A 490, 1 (1988). 

D. R. Tilley, CM. Cheves, J.L. Godwin, G.M. Hale, 



H.M. Hofmann, J.H. Kelley, C.G. Sheu, and H.R. Weller, 
Nucl. Phys. A 708, 3 (2002). 

M. Pervin, S. C. Pieper, and R. B. Wiringa, Phys. Rev. 
C 76, 064319, (2007). 

S. C. Pieper, R. B. Wiringa, J. Carlson, Phys. Rev. C 
70, 054325 (2004) 

S. C. Pieper, V. R. Pandharipande, R. B. Wiringa and 
J. Carlson, Phys. Rev. C 64, 014001 (2001). 

C. Forssen, E. Caurier, P. Navratil, Phys. Rev. C 79, 
021303 (2009). 

P. Navratil, J. P. Vary, W. E. Ormand and B. R. Barrett, 

Phys. Rev. Lett. 87, 172502 (2001) 

U. van Kolck, Phys. Rev. C 49, 2932 (1994). 

E. Epelbaum, A. Nogga, W. Gloeckle, H. Kamada, 

U. G. Meissner and H. Witala, Phys. Rev. C 66, 064001 

(2002) 

D. R. Entem and R. Machleidt, Phys. Rev. C 68, 041001 
(2003) 

A. Nogga, P. Navratil, B. R. Barrett and J. P. Vary, Phys. 
Rev. C 73, 064002 (2006) 

D.R. Tilley, J.H. Kelley, J.L. Godwin, D.J. Millener, 
J.E. Purcell, C.G. Sheu, and H.R. Weller, Nucl. Phys. 
A 745, 155 (2004). 

G. Audi, O. Bersillon, J. Blachot, A.H. Wapstra Nucl. 
Phys. A 729, 3 (2003). 

R. Roth, P. Navratil, Phys. Rev. Lett. 99, 092501 (2007). 

R. Roth, Phys. Rev. C 79, 064324 (2009) 

T. Dytrych, K. D. Sviratcheva, C. Bahri, J. P. Draayer, 

J. P. Vary, J. Phys. G 35, 095101 (2008). 

T. Abe, P. Maris, T. Otsuka, N. Shimizu, Y. Utsuno and 

J. P. Vary, in New Faces of Atomic Nuclei, AIP Conf 

Proc 1355, 173 (20111: larXiv:1107.1784l 

T. Abe, P. Maris, T. Otsuka, N. Shimizu, Y. Utsuno and 

J. P. Varv. larXiv:1204.1755l 

Wolfram Research, Inc., Mathematica, Version 8.0, 
Champaign, IL (2010). 

G. R. Satchler, Direct Nuclear Reactions, Clarendon 
Press, Oxford (1983). 



